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We present a detailed analysis of the features and capabilities of Wakefield-Induced Ionization (WII) injection 
in the blowout regime of beam driven plasma accelerators. This mechanism exploits the electric wakefields to 
ionize electrons from a dopant gas and trap them in a well-defined region of the accelerating and focusing wake 
phase, leading to the formation of high-quality witness-bunches 1 . The electron-beam drivers must feature 
high-peak currents (7° > 8.5 kA) and a duration comparable to the plasma wavelength to excite plasma waves 
in the blowout regime and enable WII injection. In this regime, the disparity of the magnitude of the electric 
field in the driver region and the electric held in the rear of the ion cavity allows for the selective ionization 
and subsequent trapping from a narrow phase interval. The witness bunches generated in this manner feature 
a short duration and small values of the normalized transverse emittance ( k p a z ~ k p e n ~ 0.1). In addition, 
we show that the amount of injected charge can be adjusted by tuning the concentration of the dopant gas 
species, which allows for controlled beam loading and leads to a reduction of the total energy spread of 
the witness beams. Electron bunches, produced in this way, fulfil the requirements to drive blowout regime 
plasma wakes at a higher density and to trigger WII injection in a second stage. This suggests a promising 
new concept of self-similar staging of WII injection in steps with increasing plasma density, giving rise to 
the potential of producing electron beams with unprecedented energy and brilliance from plasma-wakefield 
accelerators. 

PACS numbers: 52.40.Mj, 52.65.Rr, 52.25.Jm, 52.59.-f 


I. INTRODUCTION 

Beam driven plasma wakefield accelerators (PWFA) 2 ’ 3 
can generate and sustain accelerating gradients in ex¬ 
cess of ~ 10 GV/m over meter-scale distances. This was 
proven in a landmark experiment at SLAC 4 , where the 
energy of the electrons in the tail of a 42 GeV electron 
beam was doubled within a distance of less than a me¬ 
ter, after being accelerated in the plasma wake, which 
was excited by the head of the beam. Harnessing such 
extreme fields for the production of multi-GeV energy, 
high-brightness electron beams will enable a new genera¬ 
tion of accelerators, capable of compactly driving appli¬ 
cations e.g. in particle physics, optics, medicine or mate¬ 
rials science. 

Improved control over the injection of beams into a 
suitable phase of the plasma wake is a necessary step 
towards the production of beams from plasma-based ac¬ 
celerators with a quality comparable to those generated 
in current state of the art particle accelerators. The first 
peaked electron beam spectra were obtained from laser 
driven plasma wakefield accelerators (LWFA) 5 by uncon¬ 
trolled nonlinear wavebreaking 6-8 . Subsequently, the im¬ 
plementation 9-15 of earlier proposed 16-19 controlled in¬ 
jection methods paved the way towards electron beams 
with improved tunability, stability and quality. The ex¬ 
perimental realization of injection techniques in PWFA 
is not yet as developed as in LWFA, and only recently 
PWFA has been advanced from single-beam electron 
acceleration 4,20 to the acceleration of distinct witness 
beams 21 . To obtain high-quality beams from PWFA, dif¬ 
ferent injection techniques, based on density tailoring 22 , 
magnetic-induced trapping 23 or external witness bunch 


generation 21 have been proposed. A promising approach 
is the injection of electrons in plasma by means of field- 
induced ionization of a dopant gas with appropriate ion¬ 
ization potential 24,25 . Injection from ionization, induced 
in a dopand gas by the radial electric field of the driving 
beam was first observed in an experiment 24 at FACET 26 . 
That method is highly sensitive to details of the initial 
charge density distribution and the transverse emittance 
of the driving beam and is therefore not easily controlled 
and contingent on fluctuations of the beam structure. In 
addition, the injection of excessive charge can severely 
beam load the wake and may result in the production of 
a continuous current that distorts the wakefield 2 ' instead 
of a well defined witness beam. Laser-triggered ionization 
injection 28,29 offers the possibility to control the volume 
of injection in a precise way. These methods rely on a 
high degree of spatial (~ 1 pm) and temporal (~ 10 fs) 
alignment of the laser pulses with respect to the electron 
beam in order to induce ionization in the desired phase 
of the plasma wake and to produce high quality beams. 

A novel ionization-based injection strategy, Wakefield- 
Induced Ionization (WII) injection, was recently pro¬ 
posed 1 as a method to achieve a high degree of control 
over the injection of electrons into the appropriate phase 
of the plasma wave. In contrast to the aforementioned 
methods, WII injection is neither sensitive to details of 
the beam structure, nor does it rely on additional de¬ 
vices, such as lasers, for the injection. Instead, it ex¬ 
ploits the difference of the absolute electric-field strength 
between the accelerating and decelerating regions of the 
first wakefield bucket in the blowout regime 30-32 to se¬ 
lectively ionize a small volume of a background dopant 
gas near the phase of maximum acceleration. In this 
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way the production of high-quality, ultra-short (~ fs), 
low-emittance (~ /rm), multi-GeV-energy electron beams 
from a comparably simple experimental setup is possible. 

In this work, we present a detailed analysis of the fea¬ 
tures and capabilities of WII injection. Since the method 
relies on properties which are inherent in the structure of 
beam-driven blowout plasma wakes, we first review the 
phenomenology of PWFA in the blowout regime (Sec. II). 
The most important scalings are discussed there and 
compared to data from 3D simulations with the Particle- 
In-Cell (PIC) code OSIRIS 33-35 . This analysis allows 
for a thorough determination of the requirements for 
WII injection. The injection principle is demonstrated 
in exemplary 3D simulations with OSIRIS, considering 
two realistic experimental scenarios for PWFA, namely 
the FLASHForward project at DESY (Sec. Ill) and the 
FACET experiment at SLAC (Sec. IV). WII injection 
also allows controlled beam-loading to reduce the total 
energy spread of the witness bunch (Sec. V). This results 
in the generation of compact, high-energy, low-emittance 
and high-current witness bunches, which fulfil all require¬ 
ments to again trigger WII injection as driver beams in 
substantially higher density plasmas. This brings up a 
new concept of self-similar staging, which has the po¬ 
tential of producing electron beams with unprecedented 
energy and quality in PWFA (Sec. VI). 


II. THEORY 

PWFA use relativistic charged particle beams that are 
sufficiently dense to significantly displace the plasma elec¬ 
trons from their ions by means of the Coulomb force of its 
charge. As the beam passes, the displaced plasma elec¬ 
trons are attracted back by the excess of positive ions 
left behind, and oscillate around their equilibrium po¬ 
sition, generating in this way a plasma density pertur¬ 
bation that is propagating at the velocity of the beam 
(vb ~ c). For small displacements, these oscillations are 
harmonic at a frequency given by u> p = \Jnoe 2 /meo and 
a wavelength of X p = 2’njk p = 2 ’kc/uj p , where no is the 
plasma particle density, eo is the vacuum permittivity, c 
the speed of light, and m and e are the electron mass and 
charge, respectively. High-current (/& > 1 kA), resonant- 
length (Lb ~ X p ) and narrow ( k p cr r < 1) electron drivers 
blow out all the plasma electrons from their propaga¬ 
tion path, creating an ion cavity with no plasma elec¬ 
trons inside. As can be seen in PIC simulations (Fig¬ 
ure 1 (a)), this ion-cavity (or blowout) is clearly delim¬ 
ited by a sheath of plasma electrons, which accumulate 
at a distance r; from the propagation axis. The blowout 
regime features ideal properties for the acceleration and 
transport of compact electron bunches 30 . Given the ab¬ 
sence of radial currents inside the blowout (and assum¬ 
ing cylindrical symmetry), it can be directly shown from 
Maxwell equations that inside the ion cavity the trans¬ 
verse wakefields, W± = E r — cB depend linearly on 
the radius Wj_/E 0 = k p r/2 7 and are constant along the 



Figure 1. An OSIRIS 3D simulation of a high-current 
(I b = 10 kA), moderately wide (( k p a r = 0.8)), axially sym¬ 
metric Gaussian electron beam, going through a plasma at 
the (linear) resonant length (k p a z = \/2). (a) Spatial par¬ 
ticle density, (b) Longitudinal electric field, (c) Transverse 
wakefield. (d) Electric field magnitude, (e) Wake potential. 
Red solid lines are the corresponding quantity along the on- 
axis region, except for (c) and (d), where the profile is taken 
0.1 kp 1 off-axis. The outer and inner circles represent the 
blowout radius estimations through Eq. (6b) and Eq. (9), re¬ 
spectively. The dark dotted lines in (b) indicate the model 
estimations for the maximum decelerating field in the beam 
region (right horizontal line, Eq. (6a)), the maximum acceler¬ 
ating field (left horizontal line, Eq. (8a)), and the longitudi¬ 
nal field slope around the center of the cavity (diagonal line, 
Eq. (7)). 


co-moving variable d(W±_ = 0, where the co-moving vari¬ 
able is defined as ( = z — ct (see Figure 1 (c)); while the 
longitudinal wakefield, W z = E z , is constant along the 
radius in the blowout, d r E z = 0 (see Figure 1 (b)). Here 
Eq = k p mc 2 /e is the cold non-relativistic wavebreaking 
field 36 . The longitudinal structure of E z , on the other 
hand, exhibits a substantial difference in magnitude be¬ 
tween the decelerating part, at the region of the driver, 
E b z = E z (C,b)i an d the accelerating part, at the rear of 
the blowout, E™ = E,(( w ) (Figure 1 (b)). The ratio 
of these two magnitudes is called the transformer ratio 
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R w = \E™/E b z \. The energy gain of a witness electron 
beam, placed at the rear position of the wake is thus 
given by me 2 = R w jb me 2 after energy depletion 
of an electron driver beam with initial mean Lorentz fac¬ 
tor of 7 b. Longitudinally symmetric drivers cannot excite 
plasma waves with transformer ratios R w > 2 in the lin¬ 
ear regime 5 '. However, in the blowout regime, symmetric 
Gaussian profiles can generate transformer ratios R w > 3 
(cf. Figure 1(b)), and triangularly ramped profiles 38 can 
reach even R w > 5. The Wakefield-Induced Ionization 
(WII) injection method exploits the fact that wakefields 
in the blowout regime in PWFA have transformer ratios 
significantly greater than one to induce ionization and 
trapping of high-quality electron bunches into the ex¬ 
treme accelerating fields of the plasma wake. These elec¬ 
trons originate from an atomic species with high ioniza¬ 
tion threshold (HIT) which is doped into the background 
plasma in a short axial region of the plasma target. 

The ionization process caused by static (or slowly vary¬ 
ing compared to the ionization process) electric fields of 
a magnitude sufficient to significantly deform the atomic 
potential barrier can be described by a tunnelling prob¬ 
ability 39 , which has been determined for a number of 
atomic species 40 . In Figure 2 the ionization probability 
rates (Wi) are depicted as a function of the total electric 
field E = \E\, for some different types of gases. These 
probability rates grow exponentially as soon as E ap¬ 
proaches an ionization threshold, which is defined here 
as the field £) on for which the ionization probability rate 
becomes W\ = 0.1 fs _1 . According to this definition, 
the ionization threshold is written in units of GV/m in 
Figure 2 for the set of selected atomic species. 


'(/) 
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Figure 2. Ionization probability rates according to the theory 
by Ammosov, Delone, and Krainov (ADK) 40 . The ionization 
thresholds are indicated by the dots and their numerical val¬ 
ues written in units of GV/m. The flags on the x-axis are 
placed at field values corresponding to Eq for four different 
plasma densities. 


Generally, the field configuration of the beam driver 
and wakefield in the blowout regime enables both ioniza¬ 
tion and trapping from two regions within the first wave 
bucket 25 . One region is located at the driver beam po¬ 
sition, where the combination of the radial electric field 


of the beam driver and the electric wakefield can induce 
ionization, if sufficiently high. In the following the max¬ 
imum electric field magnitude in the region of a high 
current, tightly focused beam driver is computed to es¬ 
timate the conditions for the onset of ionization in this 
region. The radial electric field in the blowout can be cal¬ 
culated via Gauss’ law, for a given charge density of the 
drive beam and with the homogeneous ion charge density 
in the blowout. For axially symmetric Gaussian beams 
(rif, = n° e~ < -’ / 2a ~~ r / 2o >), the radial electric field is then 
given by 
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where Ab{0 = fc 2 J 0 °° ^rdr is the normalized charge per 
unit length of the beam. Ab is related to the electric 
current of the beam by Ab = (fc 2 /27recno) h = 2 R/Ia , 
where I a = income 3 /e ~ 17.05 kA is the Alfven current. 
The above defined beam thus has the normalized current 
profile A b = (n°/n 0 ) ( k p a r ) 2 e - ^/ 2 ^. The first term in 
Eq. (1) corresponds to the contribution from the beam 
charge density, while the second term represents the con¬ 
tribution from the uniform charge density of ions plus 
the slope of the longitudinal wakefield inside the cavity. 
In the drive beam region, k~ 1 d^E z /E 0 -C 1 and this con¬ 
tribution to the second term in Eq. (1) can be neglected. 
Furthermore, for high-current (Ab > 1), tightly focused 
( k p a r 1$ 1) beams, the first term in this equation domi¬ 
nates over the second term and the maximum of E r can 
be obtained analytically from Eq. (1) 
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This maximum value is reached at distance r ~ 1.585 o> 
from the center of the beam’s propagation axis. Given 
this inverse proportionality of E™ ax on a ri the magni¬ 
tude of the electric field in this region is highly sensitive 
to the betatron oscillations of the beam envelope a r in 
the focusing ion-cavity 41,42 and fluctuations in the slice- 
emittance and current profile. 

The second region from which ionization and subse¬ 
quent trapping can occur is located at the rear of the 
ion-cavity, where only the wakefields induce ionization. 
While the electric field magnitude in the beam region 
depends on details of the beam’s charge distribution, the 
electric wakefield behind the driver, in the back of the 
cavity, does not depend on these details, but only on the 
geometrical structure of the blowout. An important fea¬ 
ture of the blowout regime may be noted in this regard: 
The structure of the plasma wake is mainly sensitive to 
the integrated charge per unit length of the driver, and 
not to its transverse distribution, as long as the enve¬ 
lope of the drive beam is sufficiently smaller than the 
blowout radius ay -C r/. This can be seen in Eq. (1), 
where the impact of the electric field on the plasma elec¬ 
trons in the sheath for r = r; oy is dominated by the 
normalized current profile of the drive beam Ab. This 
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implies that the wakefields are not sensitive to betatron 
oscillations and fluctuations of the driver beam distribu¬ 
tion, and remain stable. The WII injection exploits this 
fact to induce ionization by means of the wakefields only, 
while avoiding any contribution from the electric fields in 
the driver’s region. For this reason, the dopant species 
has to be chosen to have an ionization field threshold 
above the maximum value in the driver’s region E b , and 
smaller than the electric field magnitude at the back of 
the cavity E w . For high-current, moderately wide beams, 
i.e. k p a r ~ 1, and ay < rj, the radial electric beam field is 
significantly smaller than the longitudinal magnitude of 
the wakefield in regions near the propagation axis. Given 
the absence of any current, also the magnitude of the 
electric field in the back of the cavity behind the drive 
beam can be approximated by the longitudinal compo¬ 
nent only. The condition for ionization to occur only in 
the rear part of the cavity might hence be written as 

\E Z \ > E ion > \E b z \ , (3) 

with E- Ion the ionization threshold of the given HIT 
species. In order to complete the injection, ionization 
must be followed by trapping of the released electrons. 
Electrons are only trapped if they do not escape the fo¬ 
cussing channel of the blowout while gaining a sufficiently 
large velocity to co-propagate with the plasma wake at 
Vb ~ c. Considering that the electrons are released at 
rest in a phase position Q. the velocity of the ultrarela- 
tivistic wake can be only acquired if the electrons end up 
in a phase position that satisfies 

= ipi - xj) f = 1, (4) 

where ip is the normalized wake potential, related to 
the wakefields by E z /E 0 = —d<^ip/k p and Wj_/E 0 = 
—d r ip/k p . Eq. (4) is a necessary but not sufficient condi¬ 
tion for trapping, which is obtained from the constant of 
motion of electrons 24 under the quasi-static approxima¬ 
tion 45 . Thus, in order to allow trapping from ionization, 
the plasma wake must provide a minimum difference in 
normalized wake potential of A ip = 1. Figure 1 (e) shows 
the structure of ip in the central slice of the example 
simulation. The maximum potential ip m is reached at 
the center of the blowout cavity (Cm), where E z changes 
from positive to negative and the cavity radius is max¬ 
imal (r; = r m ). The minimum ip m \ n is located at the 
very end of the cavity, Cc, where the plasma electrons 
cross the propagation axis and E z changes from negative 
to positive values (cf. Figure 1(e)). The maximum po¬ 
tential difference of ip is therefore given by the integral 
A'0 max = —k p f^ c (E z /E 0 ) dC- The phenomenological 
models for the blowout regime 31,32 connect E z with the 
dynamics of the plasma electrons in the plasma sheath. A 
particularly clear differential equation for E z is provided 


d c E z 

k p E[j 
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Eq. (5) explicitly expresses the dependence of d^E z on 
the cavity radius rj, the magnitude of E z itself, and fi¬ 
nally, the longitudinal velocity pi z of the plasma electrons 
in the screening layer, and it is used for estimations in 
the context of this paper. For sufficiently high-current 
Gaussian beams ending before the cavity center, this 
model (Eq. (5)) connects both the maximum decelerating 
field in the region of the drive beam, E b = E z ((b), and 
the maximum blowout radius, r m , with the peak current 
of the drive beam, A°, 

E b /E 0 *^A°/2, (6a) 

( k p r m ) 2 ~ \J 32ttA£ ( k p a z ). (6b) 

In addition, it can be seen that, in case of a large blowout 
radius k p r m ~3> 1, the electrons in the plasma sheath 
acquire high speed in backwards direction in respect to 
the driver beam propagation. If the beam has a negligi¬ 
ble current at the region of maximum radius within the 
blowout, one finds from Eq. (5) that E z is proportional 
to ( in a wide region from the cavity center towards the 
back. The slope of the longitudinal field is defined en¬ 
tirely by the plasma sheath velocity 

d c E z _ fa z _ 1 

k p E 0 1 - p u ~ 2 ’ U 

where the last approximation is obtained in the limit of 
Pi z —> — 1. Under the assumption that E z continues with 
this linear slope from up to the end of the cavity ( c , 
it is possible to find a simple estimate for the magnitude 
of and the maximum difference in wake potential 

K/Eo « ^ , (8a) 

A^ max «^) 2 , (8b) 

where l m = <p m — C, c is defined as the longitudinal semi¬ 
axis of the ion-cavity (Figure 1 (a)). 

In this work, we compare these approximate scalings 
with results from 3D simulations obtained with the PIC 
code OSIRIS. The series of simulations consists on axi¬ 
ally symmetric Gaussian beams of fixed geometrical di¬ 
mensions ( k p a z = \/2 and k p a r = 0.1), fixed normal¬ 
ized emittance ( k p e n = 0.31) and variable peak current 
I 3 , ranging from 1 kA up to 60 kA. Figure 3(a) shows 
r m (black line) and l m (grey line) as a function of the 
peak current of the beam. The top dashed line is the 
blowout radius obtained from Eq. (6b) for k p a z = \/2 : 
i.e. 8 a/7tA°. It is worth noting that, despite the fact that 
this radius overestimates r m for relatively low currents, 
it reproduces the scaling of l rn . The second dotted line 
is defined by 
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Figure 3. Scaling of key parameters in the blowout regime 
of PWFA as a function of the peak current of Gaussian drive 
beams, (a) Transverse (r m ) and longitudinal ( l m ) semi-axis 
of the ion cavity plus the semi-axis of the trapping zone (It). 

(b) Derivative of E z with respect to £ at the cavity center. 

(c) Magnitude of the maximum E z at the back end of the 
cavity (grey line), inside of the focusing region (black line), 
in the trapping zone (blue line) and in the driver region (red 
line), (d) Maximum wake potential difference in the blowout 
cavity (grey line) and in the focusing region (black line), (e) 
Maximum transformer ratio in the cavity (grey line) and in 
the focusing region (black line). The dotted lines represent 
certain model estimations that are explained on the text. 


which is a well known scaling of the transverse blowout 
radius as a function of the driver current, for moder¬ 
ate length Gaussian drivers (k p a z ss v2), obtained from 
PIC simulations 32,41 . For high peak currents, r m ob¬ 
tained from Eq. (9) as well as r m obtained from Eq. (6b) 
asymptotically converge towards l m in the PIC simu¬ 
lations. The sheath of plasma electrons hence aproxi- 
mately follows a spherical shape behind the driver beam 
for high current drivers. The scaling for the blowout 
radius with the beam peak current in Eq. (6b) in this 
comparison is identified as an accurate estimation for 
l m (cf. Figure 3(a)). Hence, the maximum ip difference 
may be connected with the current of the driver beam 
A^max ~ \J^ A 3 , by means of this scaling. In contrast, 
the scaling for r m from Eq. (9) was used in a previous 
work 25 to express the maximum wake potential difference 
as a function of the driver’s current A ip max « A 3 , and 


hence, estimate the minimum peak current for Gaussian 
beams to allow for trapping from ionization: A 3 > 1. 
In Figure 3(d), we show the scaling of Ai/) max , as ob¬ 
tained from full 3D PIC simulations (grey line), together 
with the two estimations above. The dotted line shows 
^/47rA 3 , while the dash-dotted line depicts A 3 . Although 
slightly above, the yj A 3 scaling seems to give a better 
description of the scaling of A ip max from PIC simulations. 

To be trapped and transported, electrons not only need 
to reach the velocity of the wake, but also must remain in 
the focusing region of the plasma wave. PIC simulations 
show that the position where the cavity closes and where 
ip is minimum falls in a defocusing region, so that elec¬ 
trons cannot be transported there. The trapped orbits 
have to lie completely in the focussing region, ahead of 
the co-moving position Cf at which the field configuration 
changes from focusing to defocusing (Figure 1(c)). The 
black line in Figure 3(d) shows the maximum potential 
difference that can be reached inside the focusing region, 
i.e. ipm — ip F , which is clearly smaller than the absolute 
A'i/’max (grey line). The maximum accelerating field in 
the focusing region is given by E f = E z (C f) (black line 
in Figure 3(c)), to be compared with the absolute min¬ 
imum at the very end of the cavity (grey line). To be 
trapped exactly at Cf, an electron must be ionized at a 
position Ci > Cf, such that ipj — ipp = 1. The set of 
points fulfilling this condition form a surface that delim¬ 
its the trappable zone. Electrons released (at rest) inside 
this zone can potentially be trapped in the focusing re¬ 
gion of the blowout cavity. The blue line in Figure 3(a) 
shows It = Cm — Ci (th e trapping distance) as a function 
of 4°, as obtained from the set of PIC simulations. We 
see that a trappable zone (It > 0) only exists for suffi¬ 
ciently high peak currents 7° > 5 kA. This defines a 
lower limit for the peak current of Gaussian drivers to 
achieve trapping of electrons based on ionization injec¬ 
tion techniques in the quasi-static picture. Nevertheless, 
we would like to point out here that this limit can be 
lowered if the phase velocity of the wake is reduced when 
performing the ionization in a plasma density transition. 
In order to induce ionization in the trappable region by 
means of the wakefields at the back of the blowout, the 
magnitude of the accelerating field at Ci, Ep = E z ((i) 
(blue line in Figure 3(c)) has to be greater than the ion¬ 
ization threshold of the dopant species \E{\ > E ion . In 
addition, it is important to avoid any spurious contri¬ 
bution from the evolving electric fields in the driver’s 
region (cf. Figure 4(b)). As we saw from Eq. (1), the ra¬ 
dial electric field is sensitive to the envelope betatron 
oscillations of the beam, and it increases significantly 
when the beam envelope is reduced. Electrons ionized 
by the radial fields in the driver beam region gain trans¬ 
verse momentum while they are still slow (<C c), and 
are likely to escape the relativistic focusing cavity in the 
transverse direction. Another reason to avoid ioniza¬ 
tion in the beam region is the fact that narrow beams 
can ionize entirely the dopant species around the axis, 
leaving no atoms or molecules which can be ionized in 
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Figure 4. Evolution of the longitudinal electric field (a), 
the total electric held (b) and the ionization rate from He 
(c) in the on-axis region, for an unmatched 10 kA Gaussian 
electron driver. The dotted lines delimit the on-axis phase 
range with ip — ipp > 1. The dotted-dashed lines represent 
the end of the region where the radial Helds of the beam are 
highly defocusing for slow electrons ( E r <C 0). 


the rear wakefields. For these reasons, beams in WII 
injection should not trigger ionization themselves in re¬ 
gions near the axis. It is interesting to note that the 
beam spot size that makes in Eq. (2) equal to E b z 

in Eq. (6a) is given by k v a r ~ 0.64 \/A^, which is al¬ 
ways much smaller than the blowout radius in Eq. (9). 
The latter condition ensures that the total electric field 
is dominated by its longitudinal component in the beams 
region near the propagation axis, while mantaining the 
condition for a stable blowout ( a r < r m ). However, un¬ 
less the beam is matched to the focusing plasma channel, 
this situation does not last, and the injection must occur 
before the beam is transversely compressed due to the 
betatron oscillations within the ion channel. After the 
injection, it is preferable that the driver beam is tightly 
focused to ensure a stable propagation that mitigates the 
head erosion 44 . This sets a condition for the normalized 
transverse emittance of the beam to be smaller than the 
matched value k p e n < y/'f /2 ( k p a r ) 2 . The spatial period 
of the envelope oscillations of the beam is given by half 
the betatron wavelength 42 of the beam-particles in the 
blowout \p = yf2r)/k p . Figure 4 shows the evolution of 
E z (a) and E = \E\ (b) in the on-axis region ( k p r < 0.1), 
for the Gaussian driver beam in the example of Figure 1. 
Due to the envelope oscillations of the beam, the total 
electric field in the driver’s region oscillates with a spatial 
periodicity of Xp/2 (Figure 4(b)). On the contrary, the 
longitudinal fields (in Figure 4(a)) are stable during the 
propagation. After the beam is transversely compressed, 
the region with a high ionization rate moves from the 
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Figure 5. Schematic of the plasma-cell setup assumed as a 
basis for the OSIRIS 3D simulations. A thin jet of a neutral 
hydrogen/helium gas mixture upstream of a Hat-top hydrogen 
gas target assures localized WII injection. The hydrogen is 
pre-ionized with a laser while the helium remains non-ionized. 


back of the cavity to the front of the beam (Figure 4(c)). 

Figure 3(c) depicts the on-axis values for the longi¬ 
tudinal electric field at the smallest possible co-moving 
position from which electrons can still be trapped E\ = 
E z (£/) (blue line) and the maximum longitudinal electric 
field in the beam region E b (red line). The longitudinal 
electric field in the beam region according to Eq. (6a) is 
represented by the red dotted line, and closely follows the 
curve from the simulation data. The grey dotted line in 
Figure 3(c) is the result of extrapolating the linear ten¬ 
dency of E z at the cavity center (Eq. (7)) up to the end 
of the cavity £ c . PIC simulations show that the difference 
in magnitude between E z and E b z grows with increasing 
peak current of the driver, hence allowing a wider range 
of ionization thresholds for which the ionization in the 
beam region is suppressed but for which the ionization 
and trapping in the blowout is still possible. The ap¬ 
propriate dopant species for WII injection must have an 
ionization threshold \E{\ > E[ on > \E b \ (blue and red 
lines in Figure 3(c)). Such an ionization threshold can 
only be obtained using high-current beams 7° > 8.5 kA. 
The criterion for the current corresponds to a normal¬ 
ized current of A 3 ~ 1, and is independent of the plasma 
density. In addition, Figure 3(b) depicts the slope of E z 
around the cavity center (grey line), which converges to 
the value suggested by Eq. (7) for high-current drivers. 
Figure 3(e) shows the overall maximum transformer ra¬ 
tio (grey line) and the maximum transformer ratio for 
witness beams placed within the focusing region (black 
line) in the PIC simulations. 

The theoretical framework for WII injection described 
above is independent of the choice for the plasma den¬ 
sity. However, the choice of a specific plasma density 
with a specific Eq requires the choice of an appropriate 
dopant species with suitable ionization threshold(s) (Fig¬ 
ure 2). Ideally, the plasma density is chosen such that the 
drive beam generates wakefields with the greatest possi¬ 
ble transformer ratio. This is the case at the so-called 
resonant length, which in the linear regime 45 is given by 
k p <J z = \[ 2 . The horizontal lines in Figure 1(d) indicate 
the ionization thresholds for H, He, and Ne in case of 
Eq = 105 GV/m, which corresponds to a plasma density 
of n o = 1.2 x 10 18 cm -3 . For this particular example, 
the first ionization thresholds of He and Ne (solid or- 
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ange and green lines in Figure 1(d), respectively) fulfil 
the condition expressed in Eq. (3) and can be used for 
WII injection. 

In order to exemplify and further analyze the charac¬ 
teristics of WII injection and the trapped bunches, results 
from PIC simulations on two experimental setups will be 
presented in the following. These are the FLASHFor- 
ward project at DESY and the FACET facility at SLAC. 
The experimental setup considered is sketched in Fig¬ 
ure 5; The plasma is pre-created by a laser pulse with 
an intensity sufficiently high to fully ionize a gas with a 
low ionization threshold (LIT), e.g. hydrogen along the 
propagation axis of the drive beam. In addition, a micro¬ 
nozzle 45 fed by the same LIT gas doped with a tunable 
concentration of a high-ionization-threshold (HIT) gas 
(e.g. helium), is positioned in the vicinity of the gas-cell 
entrance. Because the electrons to be injected in the wake 
by means of Wakefield-Induced Ionization will originate 
from the HIT dopant gas, it must remain non-ionized af¬ 
ter the passage of the ionization laser and the electron 
driver beam. This technology for the plasma cell 4, also 
allows the density profile of the gas components to be 
controlled by means of several gas inlets with tunable 
pressure. The density profile configuration for WII injec¬ 
tion can be simple and we only consider the case in which 
the density of the LIT gas is approximately uniform all 
along the main channel and the micro-nozzle. The doped 
gas jet emerging from the nozzle can spatially be con¬ 
fined to a distance which is significantly shorter than the 
betatron wavelength, so that an unmatched beam expe¬ 
riences maximum transverse compression downstream of 
the doped gas nozzle. 


III. FLASHFORWARD AT DESY 

FLASHForward (Future-oriented wakefield-accelerator 
research and development at FLASH) is a plasma- 
wakefield acceleration project at DESY, which utilizes 
high current and short electron beams from the versa¬ 
tile FLASH accelerator. The standard operational mode 
of the FLASH accelerator offers electron bunches with 
approximately Gaussian longitudinal (a z ~ 20 /im) and 
radial (ay = 10 /im) profiles, transverse normalized emit- 
tances of e n ~ 1 pm, and an energy of 1 GeV with a rela¬ 
tive spread of 0.1 %. Such bunches are to be compressed 
even further in the beam-extraction arc upstream of the 
FLASHForward plasma cell, to reach peak currents of up 
to lb ~ 10 kA in a more compact size (cr z ~ 7 pm and 
oy = 4 pm). Operating this current-enhanced FLASH¬ 
Forward driver at the resonant length ( k p a z = \/2) re¬ 
quires a plasma density of no = 1.2 x 10 18 cm -3 . This 
results in parameters for the electron driver and the den¬ 
sity that correspond to the parameters used for the PIC 
simulation in Figure 1 for a 10 kA Gaussian driver. As 
pointed out at the end of Sec. II, at this plasma density 
the first ionization levels of both He and Ne match the 
condition for WII injection (Eq. (3)) with a 10 kA beam. 
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Figure 6. OSIRIS 3D simulation of a FLASHForward-type 
electron beam (h, = 10 kA), traversing a plasma at the res¬ 
onant density (no = 1.2 x 10 18 cm -3 ), (a) Spatial particle 
density, (b) Ionization rate for He according to the ADK 
model. The contours in panel (b) show profiles of the equipo- 
tential surfaces. These surfaces are depicted in steps of 0.2 
starting from the minimum value inside the focusing region 
(i/'f). The vertical dotted line indicates the beginning of the 
He region. 

Figure 6(a) shows a snapshot of the simulation when the 
beam traverses a Lpe = 60 pm diameter gas-jet doped 
with He at npe = 0.01 no concentration. Figure 6(b) 
shows the expected ionization rate from the outer He 
electron in a configuration with a 10 kA peak current 
beam operating in a plasma of no = 1.2 x 10 18 cm -3 . 
Because the gas streams from right to left at c in this 
co-moving picture, the ionization rate is computed tak¬ 
ing into account the He electrons that have been already 
ionized 

m) = wi(o(i-m)), (io) 

where Wi(C) = Hj(E(())), the species ionization rate ac¬ 
cording to the ADK model 40 and 

Pi (0 = c- 1 [ Wi(0 dc', ( 11 ) 

the ionization probability in the co-moving frame. The 
regions with a high ionization rate will be denoted as 
the ionization volume. The volume from which injection 
is possible is thus determined by the intersection of the 
ionization volume (E > £) on ) with the volume satisfying 
the trapping criterion (ip — ipp > 1). The contours of 
ip — ipF are drawn in Figure 6(b), showing that only a 
narrow and well-defined region around the propagation 
axis satisfies the conditions for simultaneous ionization 
and trapping. Owing to the small size of this injection 
volume, the trapped orbits will also be confined in a small 
transverse region and short phase interval during the ac¬ 
celeration process. These are prerequisites for the gener¬ 
ation of intrinsically very short and low-emittance elec¬ 
tron bunches. The ionization rate increases quickly as 















soon as the wakefield exceeds the ionization threshold, 
thereby defining a narrow spike of high ionization rate 
in the on-axis region (red line in Figure 6 (b)). Electrons 
originating from this location Q are trapped at com¬ 
plying with V’(Ci) — VKC/) = 1- Because the V’ contours 
are closer to each other when the electrons approach the 
trapped positions, the longitudinal distribution during 
the acceleration of the injected electron bunch is nar¬ 
rower than the initial distribution. In the considered ex¬ 
ample, the rms of the initial ionization volume is signif¬ 
icantly smaller than the plasma skin depth k ~ 1 . Due to 
the invariance of E z along the cavity radius, the volume 
of injection extends transversely up to 77 . Furthermore, 
the radial electric field increases linearly with increasing 
radius (Eq. (1)) and thus, the ionization rate is greater 
in positions close to the transverse boundary of the cav¬ 
ity. However, electrons ionized close to the boundary 
are, either out of the trapping region or escape trans¬ 
versely due to sizeable transverse momentum gain from 
the positive transverse electric field in the most off-axis 
regions. The initial injection volume may therefore be 
considered as a thin disc centred on axis which extends 
up to a radial position r max that is always smaller than 
the cavity radius at that point (cf. Figure 6 (a)). The 
finally injected witness bunch is composed of electrons 
ionized from within this volume. Electrons belonging to 
the same final (j slice originate from different radial posi¬ 
tions along the initial ^ contour. Assuming full betatron 
decoherence for every slice, an upper estimate of the nor¬ 
malized transverse emittance can be given in terms of 
the initial transverse extent of the slice 25 e n = fc p ( 7 , 2 )/ 4 . 
Considering for simplicity, the ionization electrons to be 
uniformly distributed up to r max along the largest ipi con¬ 
tour, the estimated maximum slice emittance yields 

,max {kp r i nax) /1rA 

K P e n — 22 ' 

Typically r max is much smaller that the maximum 
blowout radius r m provided that ionization and trapping 
only happens from the core of a thin ionization slice at 
a rear position of the blowout cavity. For this reason, 
the WII injected bunches are, by construction, extremely 
compact and low emittance. As a very practical rule of 
thumb, we can write the following relation for the WII 
injected bunches 


kpd z ^ k p e n ~ 0.1. (13) 

That is, both the length and the normalized emittance 
of the injected bunch are a fraction of the plasma skin 
depth, which scales with the inverse of the squared root 
of the plasma density (k oc 1 /i/do). 

Figure 7 shows the density of the driver-witness system 
after 14.6 mm of propagation in plasma in the simulation 
discussed above, where a short (240 nm) bunch of elec¬ 
trons has been ionized and injected from the neutral He 
by means of the wakefields only. The total injected charge 
amounts to 6.43 pC, while the rms current, defined as 
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Figure 7. Results from the same PIC simulation as in Fig¬ 
ure 6, but after 14.6 mm of propagation in plasma, (a) Spatial 
particle density. Here a short bunch of 6.34 pC has been in¬ 
jected and subsequently accelerated up to more than 2.5 GeV 
energy, (b) Longitudinal electric fields. The bunch expe¬ 
riences an accelerating gradient substantially greater than 
100 GV/m. 
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Figure 8. Witness bunch properties after 14.6 mm of ac¬ 
celeration. (a) Absolute charge distribution of the bunch in 
longitudinal phase space ( p z vs. £ plane). The projection of 
this distribution in p z is depicted on the left axis, (b) Bunch- 
current profile. The relative energy spread and the transverse 
emittance are plotted for different longitudinal slices along 
the bunch. 


It ms = cQ/(V2^a z ), is equal to 3.2 kA. Note that with 
this definition, the peak and rms currents for Gaussian 
beams are identical. The injected beam has been acceler¬ 
ated over a distance of 14.6 mm, from (Q) ss —29.8 ^m, 
where E z (((f)) ss 136 GV/m (Figure 7(b)). 

Some additional properties of the trapped bunch can 
be estimated from the initial phase-space distribution 
as well. Trapped electrons with the same initial ipi 
propagate to the same final position Q, which fulfils 
7 /, — ipf = 1; they will be accelerated by the same 
field value E z j. However, each one of these slices in 
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£/ is composed of electrons ionized at different longi¬ 
tudinal positions along the He column, and therefore 
are accelerated over different times, producing a finite 
spread in longitudinal momentum in every slice given by 
A p z j ~ — eE z j Ln e , which at the average position of 
the bunch gives A p z j ~ 8 MeV/c. This is identified as 
the main contribution to the sliced energy spread of the 
beam, as other sources of thermal spread (e.g. helium 
temperature) on the initital population of injected elec¬ 
trons are considered to be well below MeV levels. Be¬ 
cause this initial sliced energy spread due to the finite 
width of the dopant jet is not expected to grow during 
the acceleration process, the sliced relative energy spread 
of the bunch = Ap z j/( r yf)mc can reach values 

well below 1% at 1 GeV electron energies. On the other 
hand, the total relative energy spread is proportional to 
the variation of E z along the bunch, which in case of 
negligible beam loading and sufficiently short bunches is 
approximately given by Ay/(7) ~ (d^E z j / E z j) cr z j. 
This quantity does not depend either on the beam’s 
electron energy or on the plasma density, but only on 
the product of the local E z slope and its length. Us¬ 
ing Eq. (7) to approximate the magnitude of the E z 
slope, and Eq. (12) for the bunch size, A~f/(j) « 5 % 
for WII injected witness bunches. From Figure 7(b), 
d(E z (((f}) S3 28.9 (GV/m) pm -1 , and Aj/^) « 8.8%. 

The properties of the simulated injected bunch after 
acceleration over a distance of 14.6 mm are summa¬ 
rized in Figure 8. The longitudinal phase space (Fig¬ 
ure 8(a)) exhibits a thin curved chirp with an average 
energy of ~ 2.55 GeV and a projected relative energy 
spread of 8.8%. The horizontal dotted lines in Figure 8(a) 
delimit the full width at half maximum (FWHM) for the 
projected energy spectrum on the left. The value of rel¬ 
ative energy spread (~ 2 %) in FWHM is given in the 
figure. Additional bunch properties can be seen in more 
detail in Figure 8(b). The current profile has a maximum 
approximately at the center of the bunch of ~ 5 kA. 
The relative energy spread (< 0.3%), and the normal¬ 
ized transverse emittance (< 300 nm) are shown for dif¬ 
ferent slices in (. The slice width (0.005 k” 1 ) is chosen 
short enough such that a further reduction of this width 
would not significantly change the values of the sliced 
properties, in order to appropriately reflect their uncor¬ 
related values. The maximum uncorrelated normalized 
emittance along the bunch is fc p e™ ax = 0.06, which may 
be connected with the injection volume maximal radius 
as discussed above. Knowing the emittance in the simu¬ 
lation, the maximum ionization radius can be estimated 
by using Eq. (12), r max ~ 0.86 k~ x = 4.2 pm, which 
matches well with the simulations (cf. Figure 6(a)). 


IV. FACET AT SLAC 

The next example considers electron bunches similar to 
those provided by FACET at SLAC. These beams 26 can 
be approximated by Gaussian longitudinal (a z = 14 pm) 
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Figure 9. OSIRIS 3D simulation of a FACET electron beam 
(lb = 23 kA), traversing a plasma with no = 0.5 x 10 18 cm -8 
and resonantly exciting a plasma wave, (a) Spatial parti¬ 
cle density, (b) Ionization rate map for He according with 
the ADK model. The contours in panel (b) outline the iso¬ 
potential regions. They are drawn in steps of 0.2 starting from 
the minimum value inside the focusing region (i!>f)- The ver¬ 
tical dotted line indicates the beginning of the He region. 


and transverse (cr x = a y = 10 pm) profiles with peak 
currents of 23 kA, transverse normalized emittances of 
£„,! = 50 pm and e n ^ y = 5 pm, and an energy of 23 GeV 
with a relative spread of 1%. The most important differ¬ 
ence with respect to the FLASHForward case in Sec. Ill is 
the significantly higher peak current of the driver beam 
A 6 = 2.7. Figure 9 shows a simulation of the FACET 
beam transversing a no = 5 x 10 17 cm 3 density plasma 
doped with He at nn e = 0.002 no, within a gas jet of 
Lpe = 100 pm diameter. The fact that the blowout is 
bigger in this case, implies that the maximum difference 
of the wake potential Aip max is bigger in the focusing 
region of the ion cavity (cf. Figure 3). This allows for 
more flexibility for the selection of the HIT species, pro¬ 
vided that the difference between El and E b z is equiv¬ 
alently greater. This also means that the HIT species 
can be chosen in a way that ionization happens closer 
to the end of the cavity, where the blowout radius gets 
smaller so that the injection volume is more uniform and 
constrained. Figure 9(b) shows the ionization volume 
within the FACET beam blowout, on top of the wake 
potential contours. Here, the ionization volume is ba¬ 
sically a part of the trapping region, but, as discussed 
before, the electrons with large initial offsets with re¬ 
spect to the axis escape the blowout region. A sufficient 
condition for trapping of an electron is that it reaches a 
ipf contour before escaping the blowout in straight back¬ 
ward propagation. The maximum radius fulfilling this 
specific condition is r max ~ 1.6 k~ 1 = 12 pm (cf. Fig¬ 
ure 9(b)). This allows the volume of injection to be 
estimated Vjnj — ^max ACion and hence of the total 
trapped charge Q He — —eriHe T^max -^He = 7.2 pC. 
Figure 10 shows a short bunch of 0.8 pm (rms) length 
and 8.77 pC charge, propagating with an accelerating 
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tance (< 1.5 /jm) are shown for different slices in 
Considering r max = 12 /Jin and Eq. (12), the estimated 
normalized emittance is 1.4 /j m, in excellent agreement 
with the simulation. Further acceleration of the captured 
beam is possible until the driver is energy-depleted. A 
simple estimate yields the maximum achievable bunch 
energy when considering the decelerating gradients ex¬ 
perienced by the driver beam of 50 GV/rn (cf. Fig¬ 
ure 10(b)). This limits the driver propagation distance 
to ~ 46 cm, and hence, the maximum energy of the in¬ 
jected beam to about 46 GeV assuming an acceleration 
of the trailing bunch at a continuing rate of 100 GV/m. 


Figure 10. Same simulation as in Figure 9, but after a dis¬ 
tance of 20.8 mm of propagation in plasma, (a) Spatial parti¬ 
cle density. Here a short bunch of 8.77 pC has been injected 
and subsequently accelerated up to more than 2.5 GeV en¬ 
ergy. (b) Longitudinal electric fields. The bunch experiences 
an accelerating gradient largely beyond 100 GV/m. 
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Figure 11. Witness bunch properties after 20.8 mm of ac¬ 
celeration. (a) shows the absolute charge distribution of the 
bunch in longitudinal phase space ( p z vs. £ plane). The pro¬ 
jection of this distribution in p z is depicted on the left axis, 
(b) displays the bunch-current dependence on £. The relative 
energy spread and the transverse emittance are plotted for 
different longitudinal slices along the bunch. 


field of E z « 130 GV/m (Figure 10(b)). At this point 
the injected bunch has been propagating for 20.8 mm in 
the plasma wake. The bunch properties are shown in 
Figure 11 in more detail. The longitudinal phase space 
(Figure 11(a)) exhibits a linear chirp with an average en¬ 
ergy of ~ 2.6 GeV, a total relative energy spread of 6% 
and a relative energy spread in FWHM of 2.7%. The 
sliced bunch properties can be seen in more detail in 
Figure 11(b). The current profile has a maximum at 
the tail of the bunch of ~ 1.5 kA and linearly decays 
towards its front (Figure 11(b)). The relative energy 
spread (~ 0.3%), and the normalized transverse emit- 


V. BEAM LOADING 


Ionization-based injection methods have the ability of 
tuning the amount of injected charge by means of the 
dopant species concentration and/or the length of the 
injection section (approximately the gas-jet diameter in 
the here considered case). In the two previous examples 
the charge of the injected beam was chosen to be low, 
so that the beam barely deforms the wakefields, and the 
energy chirp basically is dominated by the local shape of 
the unloaded wakefield E z . In general, an electron beam 
injected in the accelerating region of the blowout cavity 
deforms the slope of the longitudinal field with respect 
to the unloaded case. For a given witness-beam current 
profile A w , this effect can be calculated from Eq. (5). 
In this discussion, we assume that the bunch is trapped 
close to, but before the end of the cavity. 

As discussed in Sec. II for high-current driver beams, 
the velocity of the electrons in the plasma sheath ap¬ 
proaches the speed of light in the backwards direction 
(Az ~ — c) when electrons in the sheath reach the point of 
maximal blowout radius. From this point on, the plasma 
electrons in the sheath are accelerated in the forward di¬ 
rection and the velocity (third) term in Eq. (5) only dom¬ 
inates at the very end of the bubble, when the electrons 
velocity approaches the speed of light in the right direc¬ 
tion (Pi z « c). In the following discussion, the sheath 
velocity term in Eq. (5) was implicitly neglected, consid¬ 
ering that the driver bunch is situated at a point where 
the electrons in the plasma sheath do not have a suffi¬ 
ciently high positive longitudinal velocity to contribute 
significantly. With this approximation Eq. (5) reads 

d c E z = 2A W 4 

k p E 0 ( k p ri ) 2 (k p ri) 2 


E z _ 

E 0 


(14) 


This equation allows the current needed to flatten the 
slope of E z along the beam to be estimated. Assuming 
perfect flattening, e.g. the left hand side of Eq. (14) to 
be identically zero gives 


A w = 2 


'El 

. Eq _ 


(15) 


where A w refers to the normalized charge per unit length 
of the witness bunch, and E™ is E z at the position of 
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the witness bunch £ w . This is a local criterion on the 
slope of E z . The shape of E z along the bunch is ob¬ 
tained through integration of Eq. (14). This has been 
done by M. Tzoufras et al. 48 for certain bunch-current 
profiles that permit for an analytical solution. Of par¬ 
ticular interest are bunch profiles that totally flatten the 
field slope along the bunch. These are found to be trape¬ 
zoidal profiles with the maximum located at the leading 
edge of the beam at C, w , linearly decaying towards the 
tail 


A 


W 
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'El ] 4 
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rpw 

fro 


(16) 


Generally, the witness bunches generated by WII injec¬ 
tion do not have this ideal trapezoidal shape for an opti¬ 
mally tailored beam loading. In the following we explored 
how the energy chirp of the witness bunch changes in re¬ 
spect to the unload case in Fig. 11, when the injected 
charge is increased. Figure 12, presents the parameters 
of injected bunches for a series of simulations identical 
other than the concentration of the dopant He in the 
gas-jet. The simulation results show witness beam pa- 
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Figure 12. Injected bunch properties after a distance of 
5.8 mm of acceleration as a function of the dopant He con¬ 
centration. (a) Total charge, (b) Current (rms). (c) Normal¬ 
ized emittance. (d) Length (rms). (e) Relative energy spread 
(FWHM). (f) Average longitudinal momentum. 


rameters at a distance of 5.8 mm downstream of the lo¬ 
cation of the injection. This is substantially before en¬ 
ergy depletion of the driving beam occurs, and the energy 
gain of the witness beam is therefore limited to about 
600 — 750 MeV. The data for the lowest He density, 
nHe = 0.002 no = 10 15 cm -3 , in Figure 12 corresponds 
to the bunch obtained in the simulation presented in 
Sec. IV. As the concentration of He is increased, both 
the total charge and the rms current of the bunch in¬ 
crease accordingly (see Figure 12(a) and (b)). In con¬ 
trast, the normalized emittance e n (Figure 12(c)) and 
the rms bunch length a z (Figure 12(d)) are not changed 
significantly for increasing bunch charge. This results 
from both the normalized emittance and the rms length 
of the injected bunch depending primarily on the volume 
of injection, so that increasing the He concentration does 
not have a severe impact. On the contrary, the projected 
energy spread of the bunch in full width at FWHM (Fig¬ 
ure 11) is dominated by the slope of E z along the bunch, 
which in turn depends on the current profile of the in¬ 
jected beam (Eq. (14)). Therefore, a clear reduction of 
the projected FWHM energy spread can be observed (see 
Figure 12(e)) for increasing charge of the injected bunch. 
The energy spread is minimal for a dopant density of 
about riHe = 0.025 no (cf. Figure 12(e)). For this case, 
the injected bunch has a total charge of 100 pC, and a 
rms current of 7 rms ss 15 kA (Figure 13(d)). Both the 
bunch length o z = 0.8 /mi and the normalized emittance 
e n = 1.3 /tm remain at the same magnitude as in the un¬ 
loaded case. If further increased, the charge of the bunch 
starts to overload and severely deform the E z slope. Fig¬ 
ure 13 shows the energy chirp for the unloaded case with 
nue/no = 0.2 % (Figure 13(a)), a partially loaded one 
with n He /uo = 1.0 % (Figure 13(b)) and the best loaded 
case with riHe/uo = 2.5 % (Figure 13(c)). More details 
on the current profile, the sliced energy spread and nor¬ 
malized emittance can be seen in Figure 13(d), for the op¬ 
timum beam loading case. The peaked projected energy 
spectrum in Figure 13(c) has a relative energy spread (in 
FWHM) of only 0.7 %, while in the unloaded case of Fig¬ 
ure 13(a), the same quantity amounts up to 4.5 %. This 
means a reduction of ~ 85% for optimal beam loading 
in respect to the unloaded case. The longitudinal phase 
space correlation is small for the high-current part of the 
beam, which is accelerated approximately at a constant 
rate along the intra-bunch axis E z « 117 GV/m. 

We compare this result with the theoretical estimations 
given in Eq. (15) and Eq. (16). Using E™ « 117 GV/m 
in Eq. (15) leads to a beam current of I w ss 50 kA, while 
using Eq. (16) with k p r m = 3.3 (Eq. (9)), one obtains 
I w ss 34 kA. The current profile of the bunch in the PIC 
simulation (Figure 13(d)) has peak currents on the order 
of 30 kA, which is not far from these analytical estima¬ 
tions. Based on these results, and on the applicability 
of the scalings for k p r m and E h z (Eqs. (9) and (6a)) with 
the peak current of the drive beam, we find the following 
expression relating the peak current of the witness and 
the driver beam with the actual driver-to-witness energy 
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Figure 13. Longitudinal phasespace of the witness bunches 
after acceleration over a distance of 5.8 mm, for the unloaded 
case (a), a partially loaded one (b) and the optimized beam 
loading case (c). The projection of these distribution in p z are 
depicted on the left axis. The bunch-current profile, the rel¬ 
ative energy spread and the transverse emittance are plotted 
in panel (d) for various longitudinal slices along the optimized 
bunch. 


transformer ratio R w 


h 



(17) 


Since WII injection requires high-current drivers with 
h > 8.5 kA and the witness bunches are injected at 
a position yielding a high transformer ratio, Eq. (17) 
implies that the produced bunches feature currents of 
tens of kA resulting in an optimal beam loading and 
the consequent reduction of the energy spread. Since 
in addition the WII injected bunches are naturally of low 
emittance (Eq. (13)), the latter result implies that the 
WII produced bunches have a high normalized brightness 
B n oc which is comparable and might be superior 

to electron beams generated in state of the art facilities 
for FEL application (LCLS, XFEL, etc.), when operating 
at plasma densities of no = 10 18 cm -3 or higher. 


VI. SELF-SIMILAR STAGING 


As a consequence of this study, provided that the WII 
injected witness bunches are short, of low-emittance and 
of high-current, they fulfil the essential requirements to 
drive strong plasma waves in a higher density plasma. 
The field configuration in this new plasma stage, in which 
the witness beam acts as the driver, may again trig¬ 
ger WII injection from an appropriate dopant species. 
Given the driver in the first plasma stage was near to 
the resonance condition, k p a z ~ 1, and the length of the 
produced witness is k p a z ~ 0.1 (Eq. (13)), the density 
in the next plasma stage, in which the witness drives 
a resonant blowout wake, needs to be about 100 times 
greater than in the first stage. We can therefore for¬ 
mulate the following relation n' 0 /no = (al/cr™) 2 ~ 100, 
where no and n ' 0 are the densities of the first and second 
stage, respectively. Owing to the increased density, the 
witness bunches generated in the second plasma stage 
(where the previous witness bunch acts now as driver) 
will have a dramatically reduced size and normalized 
emittance. Using the scaling of the witness bunch prop¬ 
erties (Eq. (13)), one can estimate a reduction of one 
order of magnitude in size and normalized emittance of 
the second witness beam with respect to the first, and 
an about two orders of magnitude increased normalized 
brightness ( B n oc / 6 /e^). In addition, due to the high 
transformer ratio, the energy per electron of the second 
witness can be double or thrice the one in the first wit¬ 
ness. This means, an increase of up to one order of mag¬ 
nitude in the energy per electron of the second witness 
with respect to the first driver. This is a novel concept of 
self-similar staging in PWFA, in which every next stage 
operates in a significantly increased plasma density. 

We have explored the validity of this new concept of 
staging using PIC simulations in the same manner as we 
did in Secs. Ill and IV. In this case we used as driver an 
electron beam with exactly the same statistical proper¬ 
ties as the witness beam in the beam-loaded FACET case 
(Figure 13 (c) and (d)). The field structure of this driver 
operating in a plasma with no = 4 x 10 19 cm -3 density 
allows for WII injection from an appropriate ionization 
level at the required ionization threshold. In this par¬ 
ticular example, the fifth ionization level of neon, with 
an ionization threshold of E- lon = 557 GV/m fulfills the 
condition expressed in Eq. (3) (cf. Figure 2). The wit¬ 
ness bunches generated in this second plasma stage fea¬ 
ture ultra-short durations (100 as (rms)), and ultra-low 
normalized emittances (60 nm), while exhibiting high- 
currents (15 kA) for a proper beam loading. Provided 
that the normalized emittance of the second witness is 
about 20 times smaller than the first witness, and that 
the peak current is essentially the same, the peak normal¬ 
ized brightness of the second witness bunch in respect to 
the first is about 400 times higher. Such values would 
largely exceed those from state of the art conventional 
linear electron accelerators. Furthermore, provided that 
the beam is injected in a phase position such that the 
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transformer ratio is around 3, its final energy after two 
stages could exceed 100 GeV per electron, if the self¬ 
similar staging process is initiated by the 23 GeV FACET 
electron beam. 

The same concept of staging can be applied to electron 
bunches produced in laser-driven plasma wakefield accel¬ 
erators (LWFA), as they could have enough current 49,50 
to drive strong plasma wakes and trigger WII injection 
in a second plasma stage of increased density (PWFA). A 
similar hybrid staging concept was proposed as a way of 
producing, from a LWFA stage, driver/witness electron 
beam pairs that can operate a PWFA stage at higher 
plasma densities than driver/witness pairs from conven¬ 
tional accelerators 51 . Here, we propose to use just one 
high-current LWFA produced electron beam that drives 
strong plasma wakes at (or near) the resonant density, 
and injects its own high-quality witness bunch by means 
of the WII injection mechanism. The requirements for 
WII injection are mainly demanding on the current pro¬ 
file of the beam, but not too strict in terms of emittance 
and energy spread. This makes electron beams gener¬ 
ated in LWFA experiments suitable candidates to drive 
the WII injection mechanism in a significantly increased 
density plasma stage, for the creation of ultra-short, low- 
emittance and high-current electron beams with rnulti- 
GeV energy in a room the size of a laser lab. 


VII. SUMMARY 

The WII injection technique requires high-current 
(lb > 8.5 kA) and moderate length (Lb ~ X p ) drive 
beams, to generate plasma waves in a strong blowout 
regime. In this regime the plasma wake is capable of trap¬ 
ping electrons from ionization. The drive beams need to 
be smaller transversely in the high-current region than 
the blowout radius (k p a r < k p r m « 2y / A/) to excite the 
plasma wave most efficiently and to transport the beam, 
but wide enough to obtain a maximum radial electric field 
below the ionization threshold in the region of the driver 
( k p a r > 0.63 \/X&). Moreover, the driver’s beta function 
needs to be greater than the matched beta in the focusing 
ion channel, so that the driver beam is transversely com¬ 
pressed after the injection to ensure a stable propagation 
that mitigates the head erosion. This sets a condition for 
the beam transverse emittance k p e n < ( k p a r ) 2 . By 

a proper choice of the HIT dopant species, the wakefield 
induces injection and trapping from a restricted phase- 
space area at the back of the first plasma bucket. The 
electron bunches injected by means of the WII injection 
method have, by construction, a characteristic size and 
a normalized emittance of the order of a fraction of the 
plasma skindepth. This means sub-micron length and 
emittance for densities close to 10 18 cm -3 . Because the 
phase of injection is located near to the end of the cavity, 
the WII injected witness beams are accelerated at (or 
close to) the optimal phase of the plasma wake, where 
the accelerating gradients are the highest and the trans¬ 


verse fields are focusing. The requirements for WII in¬ 
jection comply with those for an efficient operating mode 
of the PWFA 38 , providing the best driver-to-witness en¬ 
ergy exchange (i.e. the highest transformer ratio). This 
injection strategy results in a controlled beam loading 
to flatten the accelerating field along the witness bunch 
length, which leads to a significant reduction of the en¬ 
ergy spread. 

The required witness currents for optimized beam load¬ 
ing are on the order of tens of kA. As a result, the 
WII injection can produce high-quality electron bunches 
with short pulse length, low normalized emittance, high- 
current, high-brightness and low-energy spread, with an 
energy per electron of around three times that of the 
drive beam. This can be achieved in a relatively simple 
experimental setup. 

The witness bunches generated by the above process 
fulfil all requirements to again trigger WII injection as 
driver beams in substantially higher-density plasmas. 
This new concept of self-similar staging has the potential 
to produce electron beams with unprecedented energy 
and quality in PWFA. 
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